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ABSTRACT 
The Center for Naval Analyses Computer War Game, SEALIFT, is a Monte- 
Carlo simulation designed to help study sealift capabilities in an ASW 
environment. A mathematical model of the SEALIFT game is posited to ob- 
tain expected value results approximating those of the SEALIFT game. The 
model is cast in the Fortran terminology of SEALIFT and an effort is made 
to accurately reflect the SEALIFT flow chart Tone in its development. 


Comparisons with SEALIFT results to determine the model's reliability and 


accuracy are not made. 
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1. Introduction 

Objective. In 1963-1964 the Center for Naval Analyses, CNA, of the 
Franklin Institute constructed a rather large and complex computer war game, 
SEALIFT I [a]. This game was used as an analytical tool for the Anti Sub- 
marine Warfare Study, CYCLOPS II. Due to the large computer running time 
of SEALIFT an after the fact need for a "quick and dirty" approximation to 
actual SEALIFT results was felt, and during the summer of 1964 this author 
was privileged to work with the CNA CYCLOPS study group on that approxima- 
tion. This thesis is an outgrowth of that summer's work. The objective of 
this thesis is to provide the logic for an expected value analysis of that 
portion of the SEALIFT game directed toward calculating transport ship 
losses through submarine action. This thesis is then a "first cut" model 
to approximate a model, and the analysis is executed under simplifying 
assumptions. 

Background. Functionally, the SEALIFT game pits the submarine forces 
of one opponent against the sea transport endeavors of another opponent in 
support of a time and scope limited war overseas. The major game output 
is the number of transport ships sunk during the war and, through inference, 
the effect of these losses on the sealift supply rates. 

Structurally, SEALIFT is an event store Monte-Carlo computer simulation 
and is primarily probabilistic with essentially no kinematics or geometry 
wired in. Although certain tactical procedures were assumed during con- 
struction, they may be circumvented or modified by the input parameter set. 
This set is on the order of four hundred and contains numerous decision 
Switches as well as numerical data. It is in the scope of this parameter 
set that one finds the flexibility of the game and, unhappily, the bane of 


forthright analysis. 
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The focal point of SEALIFT is the submarine-convoy encounter and of 
the thirty-two subroutines which comprise SEALIFT, five form ee substance 
of this event. The remaining routines are directed toward making the game 
more realistic by introducing other strategies and toward controlling the 
flow of the game. Figure 1 is a graphical representation of the game. A 
general and more detailed description of the game's structure and logic is 
quoted from SEALIFT I fa in Appendix I. 

Philosophy. All submarine loss events in SEALIFT except those incurred 
in the submarine-convoy encounter are analyzed here utilizing independent 
Bernoulli trial techniques to yield simple expected value results. To 
effect this result, several event-time dependent routines were omitted from 
consideration and a few were modified to simplify the analysis. 

Preliminary investigation indicated that the submarine-convoy inter- 
action was the most sensitive in the determination of transport ship losses 
and that simple techniques were not sufficiently accurate to determine these 
losses. This encounter has been recast in the framework of a Markov chain 
where each of the sequences of events which may befall a submarine upon 
convoy passage is considered a state of the Markov chain. Primary effort 
has been directed toward accurately effecting the Markov structure, and 
minimal effort toward the accurate reflection of each state within the chain. 
Each state, i.e€., sequence of events, is described via independent Bernoulli 
trial terminology, again with approximations. 

This analysis has been made with some trepidation. It may be argued 
that detailed Markov analysis is pretentious in view of the coarse Bernoulli 
approximations; nevertheless, if the results of this expected value analysis 
approximate SEALIFT results, it is hoped this method of analysis will pro- 
vide new insight for a more detailed expected value analysis of the entire 


SEALIFT game. 
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2. Mathematical Preliminaries 


Notational usage in this thesis. The following symbols are used in 


this thesis: 
Symbol 
Pr(A) 


Pr(A/B) 


E(N) 


[| 


(qs 5) 
aga, 
A 2B 
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‘Definition 
The probability of the event A. 


The probability of the event A conditioned 
on the occurrence of the event B. 


The expected value of the random variable N. 
The greatest integer less than or equal to a. 
The brackets are used, additionally, to en- 
Close a matrix when its elements are listed 


ner. 


The matrix whose ij th element is oat 


A is defined as B, or A is identically equal 
tO Be 


A iS approximately equal to B. 


The product, Ay x A, X coo X As 


The sum, AL + A, + coo + Ae 


if and only if 





Markov chains. It is generally understood that if a physical system 
can exist in a known set of stable states, say EA: ES: ieeG \tinseve oF 
countable), and if with each state there is associated a probability, 
Pr(E.) = Psy that the system is in that state (or goes to that state), then 
the apres that the system transits the states aN oe cee a is 


the product ie a e The basic assumption here is transitional indepen- 
k=0 k 


dence, i.e., the chance for transition from a given state to another does 
not depend upon the given state or upon past states. 

Some physical systems cannot be so explained but require the dependence 
of transitional probabilities upon the history of the system. The simplest 
of such dependencies is that the next state in sequence depends upon the 
present or given state only. In this case the outcome, Eis is no longer 
associated with a fixed probability Pye but every pair (E., E.) is associ- 
ated with a conditional probability, eae which denotes the chance that 
the system goes to state = given that it was in state Ee. at the previous 
step. In general, the transitional probability Pj is a function of the 


sub-subscript in the transit E. , oe, e » and the process is called homo- 
0 n 


geneous if Pai is independent of this sub-subscript. (Only the homogeneous 
case will be considered in the expected value analysis and in what follows.) 
If, in addition, the probabilities, Ps of the initial trial i are 


known, the probability that the system transits the states E. , E. , ..., E. 


Jo oY} Jn 
ne 
is the product p. x HE pe. 2 e The conditional probability of the 
Jo. ksO) Pk k+l 
n=l 
transit, given the initial state E. , is then 7[ p.. . 
Jo k=0 2k? k+l 


The transitional probabilities can be arranged in a matrix, called the 


transition probability matrix, P 





(2.0) p 


where the standard matrix subscriptsij also stand for transition from state 
E. to state Sie Together with the initial distribution, Pas the matrix 


E 


1° ~o 2 eve & It should be 


completely defines a Markov chain with states E 

noted that the elements of each row must sum to one, and that the element 

Pai denotes the chance for transition from state E. to state E. in one step. 
The natural extension is to investigate the concept of transition 

from state E. to state ze in more than one step, say inn steps. This 

transition can occur via many paths, and hence the conditional probability, 

Ae of finding the system in state Se at the n+r th step, given that it 


was in state E. at the r th step, is the sum of the conditional probabilities 


for each possible path. In particular, 


pit) op. 
ij ij 
(2) 
Pij * Pir rj 
p 


The equation 


Ze = 
(2.1) sh) 5, pine) 


ij in 2) 


and the Chapman-Kolmogorov equation 


13 nip Pri 


(2.2) em) i ya (n) 
r 


follow by mathematical induction. 


*@i11 iam Feller, An Introduction To Probability Theory And Its Appli- 
cations (New York: John Wiley and Sons, 1957) p. 348. 
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The probabilities, eas can be arranged in a matrix of identical structure 


as that given by equation (2.0). Coupled with equation (2.1), this suggests 
matrix multiplication - the ij th element of the n step matrix is the inner 
product of the i th row of the one-step matrix with the j th column of the 


(n-1) st step matrix, i.e., 


n 1 


Pp’ = Pp xp 


Similarly, equation (2.2) implies 


Also, if P is the probability of the initial state Es the unconditional 


probability of finding the system in state E. after n steps is 


J 
(n) (n) 
(253) P = L Pi * Pi; 


A particular n step transition matrix. The expected value analysis 
requires the n step matrix, pr formed from the (one-step) transition prob- 
ability matrix, P, 


been) FS Weer) ln (ee Nyasa aie 


lt 
ae) 


ij 


BO Oo -F 90 O70) 0.0 


Y gb 0 wet 


of the previous notation, etc. If A., A 


where Pa 1? Ao» 


sen ra! 


are defined such that 
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5 air. P3 Pa Ps Pé 
then 
P = A, A, A, A, 
g Z @ 
0 61g 
Z2ggo ii 
and 
(2.5) Pomeiitae 86 (nein (A. +1)A Cesar 
a 1 i 2 1 3 1 4 
g I " "4 
g g I " 
lg d p 
and in general 
" 7 n=-1 i n=l i n=l i 
p =| AD (), AT )A, ON AL )A, OF Aaa) 
2=0 i=0 isD 
g I g g 
g " I " 
g g g I 
By recursive multiplication, it follows that 
Aaah (Geer) an Ae 
1 1 2 i a 
and therefore 
n=l i n=l n=l Pee 
ee ee ee (p, + Po) x A, +1 
i=0 teak call 


Pa Pal» 


Py Pg 





Letting 


l- (p, A py 
Fi(n) =. oe eee; F(a) 
1- (p, + p,) 


it follows that 
=I, . 
¥ va Dee teeta) 4 
: af y ae 
i=0 


Py F(n) 
and substitution in (2.4) obtains 


ola) P = (p, te p,) 


1l1- (Pp) + p,) 


= 


1- (p) + p.) 


P> fa{n) 


P, Fa(n) +] 


F fn) x 


Pz Py Ps Pg Pr Pg 
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3. The Model for the Expected Value Analysis 

In what follows, a perusal of Appendix I and of SEALIFT I | is 
assumed. 

Methodology. The heart of any computer game is the flow chart logic, 
and it was via these charts that SEALIFT was submitted to analysis. Prior 
to detailed analysis, the game was examined as a unit to cull those routines 
pertinent to the submarine from those of a bookkeeping nature and other 
routines of non-interest, including those which determined convoy losses 
to non-submarine causes. Of the thirty-two routines which comprise .the 
SEALIFT game, twenty were discarded as irrelevant to the objective of this 
thesis. The remaining routines were again scrutinized for syntax as related 
to a typical play of the game. From the observed functional and logical 
arrangement of these routines, the all important abstraction from "real 
world," SEALIFT, to model was attempted. 

Concurrent with, but not independent of, this abstraction, a variable 
which would be definitive to the expected value analysis was sought. Bear-=- 
ing in mind that the motive was to evaluate losses to the cargo fleet by 
submarine offensive, the natural and obvious choice was "the total expected 
number of torpedoes fired at cargo (vice convoy) ships during the war." As 
Stated, this number proved to be unattainable, the primary culprit being 
the time and event dependencies which permeate (are) the game. It appeared, 
however, that if the number of torpedoes a typical submarine would fire at 
cargo ships during a single cycle at sea, i.e., from home port to station 
to home port, could be determined, this number could be extended into a 
total war figure. Jhe intuitive appeal of this single cycle approach is 
that, save for the submarines initially on station, it removed, postponed 


really, the necessity for considering the element of time. It rendered 





Submarine 


\ 


Submarine Submarine 
transits transits 
barrier barrier 
\ Submarine / 
is on station 
XXXXXXXXXXXXXKXXXX 
Convoy 
( crosses subline \ 
Legend 


————»s- +: Event flow 


XXXXXX 3: Convoy-submarine interaction 


FIGURE 2 
GRAPHICAL REPRESENTATION OF THE CYCLE ABSTRACTION 


FROM SEALIFT TO THE EXPECTED VALUE MODEL 


AH 





individual submarines’ unique and time dependent lives indistinguishable 
and it collapsed the immense and complicated time-event structure into a 
fairly tractable entity. The basic cycle is represented in Figure 2. The 
model for the expected value analysis developed in this thesis is therefore 
the logic needed, using the Markov chain process, to arrive at the number 
of torpedoes a submarine would fire at cargo ships during a single cycle 

at sea. 

In developing this logic, three difficulties became immediately apparent: 

1. Aside from the three different type submarines and the four pos- 
Sible sonar boxes to which they could be assigned, the submarines initially 
on station are unique. These submarines have a probability distribution on 
their departure-from-on-station time, unlike other submarines which have a 
fixed mission time. This probability distribution necessitated a modifica- 
tion of the associated Markov chain, since prior SEALIFT results indicated 
this cycle to be a very sensitive factor. 

2. How are the results of a single cycle at sea "extended" to cover 
the total war? In particular, what assumptions are required to minimize 
the effects of a coarse approximation? How are excessively small, or large, 
convoy losses to be interpreted? 

3. How can the event and time dependencies within a single cycle be 
handled effectively? 

The third difficulty was resolved by essentially ignoring the event 
and time dependencies inherent to SEALIFT. For the most part, these depen- 
dencies proved incompatible with either the cycle abstraction or with the 
Markov chain, or both. Where possible, however, they were included within 
the framework provided by this thesis. 


The second difficulty was resolved by determining the minimum number of 
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cycles a live submarine would make during the war. The basic cycle was 
then repeated that minimum number of times with surviving submarines from 
one cycle serving as inputs to the following cycle. Rationale for this 
choice was provided by two facts. First, a priorminformation from the 
SEALIFT game indicated that where the submarine was almost immediately 
overwhelmed, the significant action took place early in the war. In this 
case, the number of cycles would prove immaterial. If, pateyew the con- 
voys were continually decimated, this expected value analysis would show 
abnormally high convoy losses, and this would actually serve the purpose 
intended by SEALIFT. 

If opposing forces were fairly balanced, successive cycles would show 
nearly identical convoy losses, again providing an indication of the progress 
of the war. In other words, under the assumptions that the omitted time and 
event dependent events are not significant, the analysis provided herein 
should reflect the trends of a corresponding SEALIFT game. This method 
of extension was also supported by the fact that, almost without exception, 
the omission of each of the dependent events effects an over-estimation of 
the number of cargo ship losses. The choice of a minimum number of cycles 
would then partially compensate for this over-estimation. 

The theoretical resolution of the first difficulty proved trivial, 
falling under the cognizance of the integral form of the theorem of total 
probability. The evaluation of this integral proved to be computationally 
messy. The functional simplicity of the probability distribution and a 
linear approximation to the step function integrand allowed a Siaaie but 
accurate approximation to this integral. Each of these difficulties is dis- 
cussed more fully in subsequent subsections. 


Having chosen a tentative model to represent the game, this model was 
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divided into its component parts and separately analyzed. The basic cycle 
could be broken into two phases, (1) the submarine-barrier penetration, 

and (2) the submarine-convoy interaction, the Markov chain of this develop- 
ment. The probabilities of the sequences of events associated with each 
phase, i.e., the barrier penetration and each state of the Markov chain, 
were analyzed via the corresponding SEALIFT subroutines or parts of sub- 
routines. The SEALIFT subroutine logic was reformulated in the structure 

of a tree or truth table. While an attempt was made to construct the trees 
to match the subroutine logic exactly, some simplication was employed. 

Basic assumptions used in the expected value analysis. Perhaps the 
most definitive assumption made in the expected value analysis is that the 
submarine-convoy encounter should be the significant event for determining 
cargo ship losses. By implication, then, the convoy losses to air attack, 
both at sea and in port, and by barrier penetration are assumed small enough 
so as not to render the submarine target-limited, vice torpedo-limited, 
and also that each convoy be initially large enough to comply with this 
restriction. Thus, Event 3 (Convoy enters barrier C1C), Event 4 (Air attack 
against convoy), Event 6 (Convoy enters barrier C2C), Event 7 (Convoy in 
delivery port), Event 8 (Convoy enters barrier C2C on return trip), and 
Event 9 (Convoy enters barrier C1C on return trip) were classified as irre- 
levant to the analysis for this reason. 

This analysis hinges, then, on the number of submarines which penetrate 
the convoys, and prior to that, the number of submarines which reach the 
subline. For this reason, all the submarine loss events were deemed impor- 
tant. In the subsequent analysis, however, it was observed that the train 
of events which could occur to one submarine was dependent upon previous 


events involving other submarines. In other words, the stochastic process 
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associated with a particular submarine was dependent upon the realization 

of the process associated with other submarines on station with it and sub- 
marines which proceeded it to station. Such dependencies are inherently 
intractable to analysis and they were ignored in the expected value model. 

By their omission, these dependencies intraduced a cumulative over-estimation 
of the number of cargo ship losses. In particular, Event 17 (Submarine enters 
variable barrier) calls the subroutine VARDET to recalculate the reduced 
barrier-detects-submarine probability in the event a submarine sinks one 

or more of the barrier ships. Following submarines would then stand a better 
chance of slipping through that barrier until the loss was replaced (Event 22 
(Replace @ barrier ship)). A concurrent effect is the torpedo usage on 
barrier ships. This expenditure directly affects the ammunition available 

to be used on cargo ships later in the submarine's cycle. This last effect 
could become more invidious during a HUK attack, Event 15 (Submarine .versus 
HUK group). In this routine the only ways a submarine could escape, if de- 
tected by the HUK Group, are if her time or ammunition runs out or if her 
tarpedoes miss their mark! In addition to introducing event dependencies, 
Pe attack event presented further theoretical difficulties in the form 
of time dependencies. Each submarine, while at sea, is subject to HUK attack 
et: HUKHR interval. This event meshes with every other event in the sub= 
marine's life and tends to hopelessly confuse expected value analysis. 

The subroutines FSEDET and ASWDET, which perform the same function for 
the forward screen escort and the ASW screen as VARDET does for’ the barrier 
ships, are called during the play of Event 19 (Submarine versus convoy). 

Like VARDET, these routines influence the train of events of following sub- 


marines. Again, parasitic torpedo usage is encountered; however, an attempt 


was made to approximate this usage in the model in the case of the forward 


13 





screen escort. 

Another, and time related, dependency was the possible monthly change 
in submarine-detects-convoy probability for the submarines stationed in the 
fourth sonar box, a portion of Event 5 (Convoy arrives at submarine line). 
This implied an incompatibility with the cycle concept. Only fortuitous 
combinations of the number of cycles of the war, duration of the war and 
constancy of the detection probability associated with the fourth sonar box 
during each cycle would permit easy inclusion of this dependency. 

Similarly, Event 13 (Strike at submarine source) was a time dependent 
event. Provision for up to ten air strikes, at times chosen through the 
input parameter set, were made in SEALIFT. These air strikes, too, are 
incompatible with the cycle concept used in the model. 

To summarize, the routines ASWDET, FSDET, VARDET, Event 15 (Submarine 
versus HUK group), Event 22 (Replace a barrier ship), Event 13 (Strike at 
submarine source), and a portion of Event 5 (Convoy arrives at the submarine 
line) were ignored by virtue of their intractability to analysis. 

Finally, Event 12 (Build submarines) provides for a monthly submarine 
biilding rate. While this event is also time dependent, it was adapted to 
the model by applying the number of submarines built during the entire war 
as an additional input to the first cycle, a striking illustration of the 
time independence of the cycle concept. 

It should be noted that, in SEALIFT, there exist at least two levels of 
logic, the internal event logic which computes losses to convoys, submarine 
losses, etc., and the logic which controls the flow of the game, i.e., com- 
putes when and how events will occur. These levels are frequently meshed 
within each subprogram, and it is with regard to the first named level that 


routines were ignored or classified as irrelevent. The second type logic 
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was used, implicitly, whenever necessary, in the analysis. In addition, 
the following were ignored as being primarily bookkeeping type routines: 
Event 1 (Build ships), Event 2 (Form convoy), Event 10 (Convoy dispersed 
at home), Event 14 (Submarine leaves port), Event 18 (Submarine arrives on 
station), Event 20 (Submarine goes home), Event 21 (Submarine arrives home), 
and Event 24 (Replace convoy ships). 

The events from SEALIFT that form the basis for the analysis are there- 
fore: Event 5 (Convoy arrives at submarine line), Event 11 (Initialize 
Submarines), Event 16 (Submarine enters constant barrier), Event 17 (Sub- 


marine enters variable barrier), and Event 19 (Submarine versus convoy \» 


3.1 Markov Chain Structure for the Expected Value Analysis 

For a given submarine, the initial Markov chain begins after that sub- 
marines first convoy has passed the subline. The initial states for this 
Markov chain are the possible sequences of events which could have occurred 
to the submarine when that convoy passed. These states may be represented 
by the arborescence, with root "Submarine on Station," exhibited in Figure 3 
where a state of the Markov chain is defined as the path from the root to a 


terminal point. Letting t. denote the states of interest, define 


aT = submarine does not detect the convoy. 
E. = submarine detects the convoy, does not penetrate the Forward 
Screen Escort (FSE), and is not killed. 

e. = submarine detects the convoy, penetrates the FSE, penetrates the 
inner screen (ASW), shoots one salvo at cargo ships, and survives 
any resulting Search and Attack Unit (SAU) action. 

> = submarine detects the convoy, penetrates the FSE, penetrates the 


ASW, shoots one salvo at cargo ships, and is killed in resulting 


SAU action. 


1l/ 





SUBMARINE IS ON STATION 


SUBMARINE DOES NOT 








SUBMARINE DETECTS CONVOY DETECT CONVOY 
E 
1 
SUBMARINE DOES NOT 
PENETRATE FSE 
SUBMARINE IS SUBMARINE IS SUBMARINE PENETRATES FSE 
KILLED BY FSE NOT KILLED BY 
FSE ; 
E, Ec. 


SUBMARINE DOES 
NOT PENETRATE 
ASW (IS KILLED) 


SUBMARINE SUBMARINE 
IS KILLED SUBMARINE IS KILLED SUBMARINE 
BY SAU ESCAPES BY SAU ESCAPES 





FIGURE 3 


ARBORESCENCE ASSOCIATED WITH THE STATES OF THE FIRST MARKOV CHAIN 
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E_ = submarine detects the convoy, penetrates the FSE, penetrates the 


ASW, shoots two salvos at cargo ships, and survives any resulting 


SAU action. 


E. = submarine detects the convoy, penetrates the FSE, penetrates the 


ASW, shoots two salvos at cargo ships, and is killed in the re- 


sulting SAU action. 


E_ = submarine detects the convoy, penetrates the FSE, and is killed 


by the ASW. 


E. = submarine detects the convoy and is killed by the FSE. 


8 


This Markov chain allows for the possibility that if a single salvo 


does not exhaust the torpedo supply of a submarine during an attack, there 


is a chance that the submarine will be unable 


to fire again at that convoy 


and will, therefore, return to station with a partial load of torpedoes to 


await the arrival of the convoys. It is assumed in this model, however, 


that two salvos will exhaust a submarine's torpedo supply. This assump- 


tion is equivalent to the assumption that the 


necessary train of events which 


would place a submarine in the position of having fired more than two single 


salvos at as many different convoys to be very improbable, and neglects this 


possibility. 


The occurrence of states EF and E imply 


that the submarine returns to 


station where she awaits the next convoy arrival, which is to say, she re- 


mains alive in this Markov chain. State E. is followed by the submarine's 


9 


attempt to return home, and is, therefore, an 


states £ EL, and E. are absorbing since 


Sara 8 


= is absorbing for another reason; it is the 
Markov chain described below. 


If the probability of the initial states 
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absorbing state. Similarly, 
the submarine is killed. State 


entry point for the second 


are denoted by PriE.) = Pas 





the one step transition probability matrix is given by equation (2.4). 

If a total of n+l convoys pass the subline during a submarine's allotted 
on-station time, then the transition matrix, P, must be raised to the n th 
power since the initial states are caused by passage of the first convoy to 
yield the n step transition matrix sha The resulting matrix, p” is given 
by equation (2.6). This matrix and the initial probability distribution, 
pin? 


Pi» suffice to determine the unconditional probabilities, » that a 


submarine will be in each of the possible states at the completion of n 
steps, (n+l) convoys having passed. The unconditional probabilities, pm, 
are obtained from equation (2.3) and from them, the probabilities of certain 
relevant events can be obtained: 


CaSO) P(A, ) 2 Pr (a submarine survives this Markov chain to attempt 


passage home/ n+l convoys have passed the subline) 


ein) pin) | a(n) 


if 2 5 
(3.1) P(A.) ¢ Pr (a submarine fires one salvo at the cargo ships in 
this Markov chain/ n+l convoys have passed the 
subline) 
(n) (n) 
= P. + PA 
(352) P(A.) : Pr (a submarine fires two salvos at the cargo ships in 
this Markov chain/ n+l convoys have passed the 
subline) 
(n) (n) 
= Pe + Pe 
(353) P(B ) = Pr (a submarine is absorbed in state E, upon passage 


of the r th convoy) 


r-1 
= aoe (p, + Po) ° PRES ee a SL 
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The second Markov chain is almost identical to the first chain and 
is concerned with those submarines which are absorbed in state E45 These 


submarines return to station having fired a single salvo of torpedoes at 


» 


the previous convoy. Assuming that a total of (n+l) convoy cross the sub- 
line, the number of steps of the second Markov chain depends upon which 


convoy caused a submarine to be absorbed in state E If a submarine is 


3° 


absorbed in state E, at the passage of the r th convoy, the (r+l)st convoy 


3 
"initializes" the second Markov chain and there are (n+l) = (r+2) + 1 = n-r 
steps in the second chain. The probability that a submarine enters the 
second Markov chain "just after" passage of the r th convoy is given by 
equation (3.3). 


The second Markov chain is described by its states: 


c 


l submarine does not detect the convoy. 


F. = submarine detects the convoy, does not penetrate the FSE, 
and is not killed. 

F_. = submarine detects the convoy, penetrates the FSE, pene- 
trates the ASW, shoots her remaining torpedoes at cargo 
ships, and survives any resulting SAU action. 

F = submarine detects the convoy, penetrates the FSE, pene- 
trates the ASW, shoots her remaining torpedoes at cargo 
ships, and ais killed in the resulting SAU action. 

F. = submarine detects the convoy, penetrates the FSE, and is 


Killed by the ASW. 


F. = submarine detects the convoy, and is killed by the FSE. 


Denoting the probabilities of the initial states of this Markov chain by 


Pr(F .) = q., the one step transition probability matrix is then 
i 


ran 





Q = 1 4) 95 95 9g 95 Ig | = 69; 
G1 92 95 94 45 4g 
: Oar 0) Wie Ceo 
UF a) peasy Say fa 
ap ia) Oto) fc) Table 


OF SOG Ue Tees] 


and the n-r step transition matrix, by analogy with equation (2.6), is 





n=-r n-r-l aL a2 Wy 42 93 44 As "6 (n-r) 
Q = (q) + q,) ; x (a5 5 ) 
q, q - GQ. 4, Ge G 
1 2 | l (q, + q,) 3) 74525 86 
Bayo ThA 


(n) 


If ips is defined as the unconditional probability that the submarine is in 


state F after n steps, then, using the theorem of total probability, 


n+l 
am) = yy Pr (a submarine is in state - after n+l convoys/ 
r= 
she entered the second chain after the r th 
convoy has passed) x Pr (she entered the second 
chain after the r th convoy). 
. (n-r) 
= y Q. <- PUBS) om = 2 
J << ~ 
r=l 
n 6 
= yy ). q. q.. ( gee » Hal 
ez) dc) 2 2 3 1 2 
where 
(o) d 
ij 7 Oj, the Kronecker delta. 
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The probabilities of the remaining events of relevance can now be ob- 


tained as: 


(3.4) P(A,) Pr (a submarine survives the submarine-convoy inter- 
action/ n+l convoys have passed the subline), 
4 * 
(n) (n) (n) (n) (n) (n) n 
= Py te Ps + Pe + RS t R + R. © p4(P, + P,) 
and 
(dea) P(AC) = Pr (a submarine fires her remaining salvo in the second 


Markov chain/ n+l convoys have passed the subline) 


(n) (n) 
R, +” RA 


In SEALIFT, the submarines may be separated into two sets, those which 
are on station at the beginning of the war and have a departure time distri- 
buted uniformly on an interval, and those submarines which arrive on station 
after the war starts and have a fixed mission time. Once the game para- 
meters have been set the arrival of each convoy is deterministic. As a 
consequence, the variable n+l is a step function of the submarines departure 
time T, i.e, n+l = g(T). Representing the equations, (3.0), (3.1), (3.2), 
(3.4) and (3.5) symbolically as 

Pr(A/n+1) = Pr (of the event A given n+l convoys have passed) 
this probability becomes Pr(A/g(T) ). Invoking the theorem of total probabil- 


ity, Pr(A) can be found as 


(3.6) Pr(A) = [ eetaetta) o Pritt 
t 


This equation is valid for both sets of submarines although trivial for 
those with a fixed mission time. For the submarine with a fixed mission time 


Pr(A) is simply the integrand of this equation evaluated at the upper limit 


2g 





of t. The evaluation of equation (3.6) is crucial in determining the ex- 
pected number of cargo ships sunk per cycle. 

This completes the construction of the Markov chain for a generalized 
submarine. To apply this structure to the SEALIFT game, the integrals Pr(A) 
and the transition probabilities Bee and she must be expressed in terms of 
of the FORTRAN variables of SEALIFT. The submarine barrier penetration must 
be analyzed in like terms. Finally, the evaluated terms must be synthesized 
to answer the basic question "What is the expected number of torpedoes fired 
at cargo ships per cycle?" and its extension "What is the expected number of 
torpedoes fired at cargo ships in the war?" 

3.2 Detailed Analysis in Terms of the FORTRAN Variables of the SEALIFT Game. 

The detailed analysis to follow will be executed in the following order, 
(1) evaluation of the integrals Pr(A), (2) evaluation of the barrier penetra- 
tion, the transition probabilities, and related expected number, (3) syn- 
thesis to determine the expected number of torpedoes fired at cargo ships 
during the war. The following list of SEALIFT I [2 FORTRAN variables is 
pertinent to the analysis, and the following convention is observed: The 
letter "I" is reserved for the type submarine index; however, it is employed 
in two ways. As a subscript, e.g, ALTIM (I), "I" can take on the values l, 
2, and 3. As a part of a variable name, e.g., NSUB'I'T, "I" is set off in 
this analysis by single quotation marks and can take on the values 1, 2 and 
N. In both usages, the values of "I" denote conventional number 1, conven- 


tional number 2, and nuclear submarines, respectively. 


Variable Name Definition 
TBCON Time between each convoy sailing (hours). 
CONSPD Speed of convoy (knots). 
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CONHR 


ONE WAY 


ULT IME 


ALTIM(1) 


DSUBTA 


SUBSPD (I) 
NSUB'I'T 
FSUB'I'OS 


el wx (ll) 


NSUB‘I' 


PBARSC(I,J) 


PBARSV (I,J) 


PVBDET (I,J) 


VAIRDET (J) 


RGDT(K,1) 


SUBLINE 


TWAR 


Time after start of war that the first convoy leaves 
home port in hours. 


One way convoy distance from home port to delivery port 
(n. miles). 


Time convoy spends in delivery port. 


Allowable time in hours for type I submarine before it 
must leave its station. 


Distance from submarine home port to submarine station 
(n. miles) e 


Speed of type 1 submarine. 
Number of type I submarines at start of war. 
Fraction of type 1 submarine on station at start of war. 


Probability that a type I submarine on station will be 
in BOX # J, J = L 000% 4. 


Number of type I submarines built per month. 


Probability that a type I submarine is sunk by the Jth 
constant barrier, J=1, oe, 4. 


Probability that a type 1 submarine is sunk by the Jth 
Variable barrier, given it was detected by this barrier, 
Ale 1 eoo, 4. 


Probability that a type I submarine is detected by the 
Jth variable barrier, with barrier at its full quota of 
ships. 


Part of detection probability of variable barrier J 

that is not supplied by ships within that barrier which 
might be sunk by a submarine. Thus, if a mixed barrier 
of aircraft and ships were used and all the ships were 
sunk, the probability of detection drived from the air 
barrier (VAIRDET(J)) would still be present. The "total" 
detection probability of a variable barrier is PVBDET(I,J) 
+ VAIRDET(J). For a ship barrier only, VAIRDET(J) should 
be set to zero. 


Maximum athwartship detection range of type 1 submarine 
against convoy type K sonar conditions 
K=l1 iff good sonar conditions 


K=2 iff poor sonar conditions 


Length of submarine line (n. miles), i.e., width of 
patrol area. 


Duration of war (hours). 
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FAIRSC 


ATSFE(1) 


PDFS(I,K) 


PDFER 


PDFCE (I,K) 


PoPSE (1) 


PREES( 15K} 


Papee(1 {Kh} 


ASALVQ(1) 


PDETU(I,K) 


VA(K) 


PBB(I,K) 


PKBB(I,K) 


Is forward escort screen made up entirely of aircraft 
(nan=sinkable)? (YES=1, NO=2). 


Number of times a type I submarine can try to repenetrate 
forward screen. 


Probability that a type I submarine in type K sonar con- 
ditions (denoted by "(1,K) submarine" henceforward) is 
detected by forward screen escort while penetrating 
screen with screen at its full quota of ships. 


Part of detection probability of forward screen that is 
not supplied by ships within that screen which might be 
sunk by a submarine. Thus, if a mixed screen of air- 
craft and ships were used and all the ships were sunk, 
the probability of detection derived from the air screen 
(PDFER) would still be present. The "total" detection 
probability of the forward screen is then PDFS(I,K) + 
PDFER. For a ship's screen only, PDFER must be set to 
zero. 


Probability that an (I,K) submarine is sunk by forward 
screen escort, given submarine has been detected by same. 


Probability that a type I submarine shoots at forward 
screen escort ship in an attempt to break contact. 


Probability that an (1,K) submarine sinks forward screen 
escort ship, given it decides to shoot at it. 


Probability that an (I,K) submarine is sunk by forward 
screen escort ship, given this submarine has already sunk 
a forward screen escort ship, i.e., flaming datum. 


Number of anti-ASW torpedoes used per salvo by a type Il 
submarine. 


Probability that an (I,K) submarine is detected by the 
ASW inner screen with direct path sonar, given it was 
not detected previously with bottom bounce sonar. 


Conversion factor reflecting improved ASW inner screen 
detection probability with direct path sonar, in type K 
sonar conditions, because of previous bottom bounce de= 
tection. 


Probability that an (I,K) submarine is detected by the 
ASW inner screen with bottom bounce sonar. 


Probability that an (I,K) submarine is killed, given 


that it is attacked or re-attacked following a detection 
with bottom bounce sonar. 
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EAE (CALS) 


PEE(1, 3) 


PABB(I ,K) 


PKBBG (I) 


Ste 17 T 


PKINS(I) 


PSCA(I,i) 


PAA(I) 


PKSC(1) 


PSUBAT (1) 


PKOSC(I) 


QTU 398 


Probability that an (I,K) submarine is attacked given 
that it has been detected with bottom bounce sonar only. 


Probability that a group from the ASW inner screen re- 
attacks a type I submarine at the i th opportunity in 
action resulting from a bottom bounce detection, 

ii 2, Se 


Probability that an (1,K) submarine attacks a group 
from the ASW inner screen in action resulting from a 
bottom bounce detection. 


Probability that a type I submarine sinks a ship in 
the group from the ASW inner screen in action resulting 
from a bottom bounce detection. 


Designator for possibility that an engagement between 
the ASW inner screen and a submarine is terminated after 
an unsuccessful submarine counter attack, e.g., due to 
quiet torpedoes used by submarine (QTC equal to zero 
implies engagement is to continue; QTC equal to one 
implies engagement is terminated. No other values are 
possible). 


Probability that a type 1 submarine is sunk by the ASW 
inner screen during action following a detection with 
direct path sonar. 


Probability that an ASW inner screen re-attacks a 

type I submarine at the i th opportunity in action re- 
sulting from a detection of sub by ship with direct 
path sonar, i=l, 2, 3. 


Probability that a type I submarine attacks the ASW 
inner screen after an unsuccessful attack on the sub- 
marine by the screen resulting from a detection with 
direct path sonar. 


Probability that a type 1 submarine kills an ASW inner 
screen ship in action resulting from detection of sub 
by ship with direct path sonar. 


Probability that a type I submarine attacks the ASW 
inner screen following an unsuccessful bottom bounce 
detection attempt by the screen. 


Probability that a type I submarine sinks an ASW inner 
screen ship without prior detection of the submarine by 
the screen. 


Designator for radiated noise level of enemy anti-screen 
ship torpedoes. (QTU equal to one implies a quiet tor- 
pedo; QTU equal to zero implies a noisy torpedo. No 
other values are possible). 
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VM(K) Conversion factor reflecting change in the ASW inner 
screen detection probability with direct path sonar in 
type K sonar conditions because of flaming datum result- 
ing from sinking of screen ship by previously undetected 


submarine. 

CSALVO(I1) Number of anti-shipping torpedoes used by type I sub- 
marine. 

PDSAU (I,K) Probability that an (1,K) submarine is detected by the 


SAV, given the submarine has sunk one or more cargo ships 
in the convoy. 


PSSAU(I,K) Probability that an (I,K) submarine is sunk by the SAU, 
given the SAU detected the submarine. 


ATTSAU(I) Does a type I submarine attack a number of the SAU group, 
given it was not sunk by this group, YES=1. 


RASW(I) Number of times a type I submarine can re-attack the 
convoy after SAU action. 


TORPS (I, ASW) Number of anti-ASW torpedoes per type I submarine. 
Note: If only one type of torpedo is used for the play 
of the game and it is desired to fire this torpedo at 
all targets, then it should [ must | be stored in TORPS 
(1,ASW). 


3.201 Evaluation of the Integrals Pr(A) 
To evaluate the integrals PrfA}, the equivalences between the first and 
second Markov chains are exploited to condense the notation. The states Eo: 


Eo: EL, and E. are seen to be equivalent to the states Fe F Fos and F 


Zs 6 


respectively, and therefore, P) = 9)» Po = Ip» Py = 46 and Pg = Pee For 
s = (p, + p>) 


and 
j-s"t 


l=-s 





f (n) 


and letting D denote the derivative operator with respect to s, the evalua-~ 


tion of the unconditioned probabilities, p(n) and ain) yield the following 


relationships: 
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p(n) 


1 aes i 

ps” = ps 

pm = p.fin), j= 3, coos 8 
Rg,” = P) P3 Ds" 

Ae = > p3 DF(n) 

as”) = q, Pp, Df(n) 

Ae = q, Py DF(n) 

oe = py Py Df(n) 

Aen = pp p, Df(n) 


Using these values, equations (3.0, (3.1), (3.2), (3.4) and (3.5) from 


section 3.1 are evaluated to yield 


(3.0) P(A.) = s™ 4p, F(n) 
(3.1) P(A) = (py +p,) fla) 
(3.2) P(A.) = (pe + Pe) f(n) 


(3.4)) P(A,) = (1 +p, D) 8" + (pe + P5q, D) F(n) 


(325) P(A.) = (q, + q,) P, D f(n) 


Generically each of the integrals Pr(A) 


Pr(A) = [ PrvAvgtt))). cd. Prt) = t) 


is a function of the integral N 


N = | 39 (t) dpe Cl = 4) 
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and the first derivative of N with respect to s, assuming the integration and 

differentiation operations can be interchanged. The distribution Pr(T = t) 

is non-trivial only for the submarine initially on station and is, in SEALIFT, 
dt 


° e ° — i, 
uniform on the interval (0, ALTIM), i.e., dPr(T=t) = ALTIM ° 


To explore the structure of the function g(t) and remain within the 
FORTRAN concept of SEALIFT, two new FORTRAN variables are defined 
TCON(T) 2 The game at which the I th outgoing convoy crosses 
the subline, 
and 


RTCON (1) The game time at which I th returning convoy crosses 


the subline. 
For those submarines initially on station, g(t) is zero for t < TCON(1), and 
steps to one at t = TCON(1). It then increases by one after each interval 
TBCON until the first returning convoy crosses the subline, at game time 
RTCON(1). Until the SEALIFT game parameters have been set, the fine structure 
of g(t) is not well defined after RTCON(1). It does, however, increase by 
two every TBCON hours after RTCON(1). Clearly then, the integrals Pr(A) 
involve extensive summations. To avoid these summations a piecewise contin- 
uous linear approximation, h(t), to the step function, g(t), is sought such 


that for g(t) = h(t), the value of N is easily found, viz., 


ai) ee f 39't) at a cat) 
z x ~ Lnr(s) dh(t) 
dt 


Incorporating the SEALIFT FORTRAN variables and the fact that in SEALIFT 
the subline is positioned half way between convoy home port and delivery port 


it can be seen that 
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_ __ONEWAY 


TCON (1) > x CONSPD 
3 ONEWAY 
RTCON (1) = 5 x CONSPD 


and further, 


TCON(T) 


RTCON (1) 


4 


TCON(1) + (I- 


CONHR 


CONHR + ULTIME 


1) TBCON 


RTCON(1) + (I-1) TBCON 


If K is defined as the number of outgoing convoys which pass the subline 


before the first returning convoy passes the subline, then K is the largest 


integer which satisfies the inequality 


TCON(K) =  RTCON(1) 


or 


(K = 1) TBCON + TCON(1) = RTCON(1) 


TBCON 


RTCON(1) = TCON(1) 


where the brackets denote "the greatest integer less than or equal to." 


Introducing a new FORTRAN variable H(I) for the expected value analysis, the 


values of H(I) and the approximation h(t) are defined as: 


H(1) ¢ h(a) = O 
d' 

H(2) = h(TCON(1) = 0) 

H(3) @ n(TCON(2) + 0) 
d 

H(4) = h(RTCON(1) = O) 

H(5) 2 h(RTCON(1) + 0) 
d 

H(6) = Hh(ALTIM) = (K 


Figure 4 depicts both a sample step 


tl 
© 


z 


1 
= 4+ TatoN 





(RTCON(1) - TCON(21)) 
= K+4 
Z 

+ +) + T5ron (ALTIM - RTCON(1)) 


function g(t), amd its linear approxi- 


mation h(t). Using h(t), the generic integral, N, becomes 
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TCON(1) RTCON(1) ALTIM 
(3.8) N = “ap + “ela + AU 
0 TCON(1) RTCON (1) 
and, by equation (3.7), becomes 
H(4) _H(3) H(6) _H(5) 
= s -s TBCON s - $s 
= TCON(1) + TBCON ints) + : Entel 


The following derivatives are of immediate interest in evaluating the inte- 


grals Pr(A), 








. Pu Titer TACO (GIN ao nD Sos Leem(ig) See N 
ALTIM ALTIM : 
Ln (s) 
SiBEORN ne i(a)ee 2 Ge es ee 
2 ALTIM ik Fo? 
and 
Ne N 
ALTIM (1 -s)] : 
(1 - s) 


In terms of these values, the integrals Pr(A) become 


(a5) PAL = Pr (a submarine initially on station survives the first 
Markov chain to attempt passage home) 
: N ALTIM - N 
~  ALTIM asp eee 
(3510) PA2 e Pr (a submarine initially on station fires one salvo at 


the cargo ships in the first Markov chain) 


a a ALTIM = N 
Te eg ed ALTIM (1 = 8s) 
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i es ee 


tf tt 
TCOWN (1) RITCON(1) ALTIM 


FIGURE 4 


A SAMPLE STEP FUNCTION, g(t), AND ITS APPROXIMATION, h(t) 
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ce 


Cae) PAI Pr (a submarine initially on station fires two salvos at 


the cargo ships in the first Markov chain) 


Sen en ALTIM = N 
Ps * Pé’  aLTIM (1 ~ 8) 


(S212) PA4 Pr (a submarine initially on station survives the inter~ 


action) 


N ALTIM ~ N 
= (1+p,D) arrm + (P5 + P3932) Artima os) 


tee 


13) PAS Pr (a submarine initially on station fires on cargo ships 


in the second Markov chain) 


Sai (accom tT —ALTIM = N 

3 4° "3 ALTIM (1 - 8) 

This completes the evaluation of the integrals Pr(A) for those submarines 
initially on station, 

The integrals Pr(A) are trivially Stieltjes integrals for those sub- 
marines which arrive on station after the war starts. The following is based 
on the assumption that the arrival time for these submarines is later than 
RTCON(1), i.€., the full convoy arrival rate of two convoys per TBCON hours 


has been achieved prior to submarine arrival at the subline. In this case, 


Pr(A) = Pr(A/ g(maximum value of t)) 
but 
ALTIM ~ submarine transit time to station 
g(max +) = | ~Aenets-Sunmarine transit time to station , 9 


TBCON ’ 


Neglecting possible delays engendered during HUK attacks and during the barrier 


: : , : d DSUBSTA 
penetration, the submarine transit time is “SUBSPD — , therefore 
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DSUBSTA 
ee SUBSPD 


Ben ainda ef 


Pr(A) = TBCON 


DSUBSTA 
SUBSPD 


TBCON 


ALTIM = 


Letting KONE : x 2 |, the following probabilities can 


be found 


ou 


(3.14) PAG Pr (a submarine not initially on station survives the first 


Markov chain to attempt passage home) 


_ KONE 1 = i 


(2215) 


(3.16) 


Cael 7) 


(3.18) 


PA? 


PA8 


PA9 


PA1O 


oO 


oOo 


tO 


ow 


Ps l-s 


Pr (a submarine not initially 
at the cargo ships in the 


ake ONE 


l-s 





(p, + P,) 


Pr (a submarine not initially 
at the cargo ships in the 


1. QKONE 


le-e-s 





(Pe + Pe) 


Pr (a submarine not initially 


marine=convoy interaction) 


(1 + py D) KONE 


Pr (a submarine not initially 


pres beyeee) oy 


on station fires one salvo 


first Markov chain) 


on station fires two salvos 


first Markov chain) 


on station survives the sub- 


1 - _ KONE 





on station fires on cargo 


ships in the second Markov chain) 


KONE 
(q Pe ge ee 
3 4 3 le-s 


This completes the evaluation of the integrals in terms of the evaluated or 
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known FORTRAN variables and the as yet unevaluated transition probabilities. 


3.2.2 The Barrier Penetration, Transition Probabilities and Related Expected 
Numbers 

Until now, it has been convenient to suppress in the analysis any nota- 
tion concerning submarine type and sonar condition. The Markov chain notation 
was difficult enough without the burden of two additional subscripts. The 
remainder of this analysis will distinguish submarines by type and by assigned 
sonar condition. It will be seen that each of the probabilities given by 
equations (3.9) through (3.18) is a function of submarine type and sonar 
condition. 

Placement of submarines. In SEALIFT, there are eight ways to classify 
submarines of a particular type by placement at the start of war: two choices 
for on station or not and four choices of sonar box. This classification is 
done by Monte-Carlo methods. Type sonar condition within each box is de- 
termined by input. For the purposes of this analysis, sonar box classifica- 
tion iS immaterial except for the restriction that the fourth sonar box is 
not included in the analysis. For the initial placement of submarines the 
expected number of type I submarines initially on station in type K sonar 
conditions is given by 

EAC) =, 9NSUBMI TT 2cPSUB ItuS x >». PY I*BOX(J) 
where the summation is over all sonar boxes assigned sonar condition K. The 
expected number of type I submarines not initially on station but assigned 
to type K sonar conditions, E3(1I,K), is derived from two sources, the fraction 
of on hand submarines not initially on station and the submarines built during 
the war. Using a nominal thirty day month, 


TWAR 


voy) XQ PY Tt BOX( J) 


E3(1,K) = (NSUB'I'T (1 - FSUB'I'OS) + NSUB'I' x 
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where the summation is over all sonar boxes assigned sonar condition K. 
Barrier penetration. A submarine will successfully transit the eight 
distinct barriers if and only if she successfully transits each barrier. In 
the SEALIFT game, the barriers are evenly divided between two types, constant 
and variable. The latter type are variable in the sense that 'the barrier- 
detects-submarine probability will vary with barrier strength. Under the 
basic assumptions used in this expected value analysis, however, each of 


the barriers is treated as a constant barrier. For constant barriers there- 


fore 
PAl1(I,J) = Pr (a type I submarine successfully transits the J th 
constant barrier) 
= 1 - Pr (she is sunk in this barrier) 
= 1 = PBARSC(I,J) 
and 
PAl2(1,J) = Pr (a type I submarine successfully transits the J th 


Variable barrier) 


1 - Pr (she is sink in this barrier) 


1 - PBARSV(I,J) x (PVBDET(I,J) + VAIRDET(J)) 
Therefore the probability that a type I submarine successfully transits the 


barriers, PO(I), is 


4 
POT )Gees 7 BPA oc ipaecr a) 
Je 


Convoy detection. In the SEALIFT game, all convoys cross the subline 
at the subline's midpoint. Detection probabilities are of the "cookie 
cutter" type where submarines are placed (at each convoy crossing) uniformly 


on the subline. The probability an (1,K) submarine detects a convoy on any 
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one trial, 1 = Pl1(I,K), is then 


2 x RGDT(K,I) 
1 -P1(I,K) = —SUBLINE 


The transition probabilities, Pp) = 94 = Pl(I,K), are thereby fully defined, 
and are dependent on submarine type and on sonar condition. 

Penetration of the forward screen eScort. [This event and the following 
events are the substance of Event 39 (Submarine versus convoy) of the SEALIFT 
game. Figure 5 is the simplified truth table (tree) of that part of Event 19 
from which the FSE associated probabilities and expectations used in this 
analysis are derived. All of these probabilities are conditioned on the 
Submarine having detected a convoy. The input parameters FAIRSC and ATSFE(1) 
are obviously not probabilities. FAIRSC is transformed into a 0,1 parameter 
in this analysis by the transformation 2 - FAIRSC, allowing its use in prob- 
ability statements. ATSFE(I) is another tactical limitation of the type I 
Submarine and this parameter is usually a "small" integer for the conventional 
submarine and a "large" integer for the nuclear submarine. If ATSFE(I) = QO, 
a type I submarine will not be allowed to attempt repenetration of the FSE 


in the event she is rebuffed. For this submarine in type K sonar conditions 


the following probabilities are obtained 


P7(I,K) = Pr (a submarine penetrates the FSE; ATSFE(I) = 0) 
= 1 — (PDFS(I,K) + PDFER) 
P9(I,K) g Pr (a submarine returns to station; ATSFE(I) = QO) 
= (PDFS(I,K) + PDFER) x (1 - PDFCE(I,K)) x (1 = (2 = FAIRSC) 
x PSFSE(T) x PREESUICK) x PSFEE(I,K)) 
and 
Pet 1) : Pr (a submarine is killed attempting FSE penetration; ATSFE(1)=Q0) 


1 - P9(I,K) - P7(I,K) 
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IS SUBMARINE DETECTED 
BY FSE? 


MES NO 
PDFS(1,K) 


+ PDFER & 


IS SUB SUNK? 


Legend: 
A: submaring penetrates FSE 


Yes NO B: submarine returns to station 
PDFCE(1,K) 


X: submarine is killed 


DOES SUB ATTACK FSE? 


YES NO 
(2 - FAIRSC) 
x PSFSE(I) 


DOES SUB SINK 
AN FSE SHIP? 


YES NO 
PFEES(I .K) 


IS SUB SUNK 


AFTER FLAMING 
DATUM? 





YES NO 
PSFCE (I,K) 


CAN SUB ATTEMPT REPENETRATION? 


NO YES 
NUMBER OF ATTEMPTS IS 
A PRE-SET DECISION GIVEN 
BY ATSFE(I) 


FIGURE 5 


ARBORESCENCE ILLUSTRATING THE SUBMARINE-FORWARD SCREEN ESCORT INTERACTION 
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In general, if ATSFE(I) =n, n attempts at penetration will be allowed, and 


P7(I1,K,n) = Pr (a submarine penetrates the FSE; ATSFE(I) = n) 
a eee ee (P9(1,K))"*> 
: 1 = P9(I,K) 
P9(1,K,n) = Pr (a submarine returns to station; ATSFE(I) =n) 
Sy yee 
and 
P8(1,K,n) = Pr (a submarine is killed attempting FSE penetration; 


ATSFE(I) = n) 


l 
; 1 - (P9(1,K)) 
Bea) eee aoe 


As derived herein these values are predicated only on the approximation that 
the FSE~detects-submarine probability has not been degraded by prior sub- 
marine attack either by the submarine under consideration (ATSFE(I) > QO) or 
by earlier submarines and that when the submarine attempts to attack the 
screen she will have sufficient torpedoes to do so. To find the expected 
number of torpedoes a submarine fires at FSE ships given she penetrates the 
screen, suppose ATSFE(I) =n, and define the following events 
= = event that a submarine penetrates the screen on the j th 
try, J = 1, ooo, nel 
D. = event that a submarine fires i salvos, each of size ASALVO(I) 
during an attempted penetration, i = 0, ..., ntl 
Then E1(1I,K) can be defined as the expected number of torpedoes fired by an 


(I,K) submarine during an attempted penetration of the FSE given she penetrates 


the ESE gamcdeATSFE(1) = n; same.) 


Sy) 





(3.19) n+l 
jee J 


1=0 P7(1I,K,n) 


CEO = ASALVOCI) 


The probability Pr (C5, D.) can be derived by the following argument: the 
event (C D.) is a sequence of j trials where the j th and last trial is 

a non-detection by the FSE. This trial is preceeded by i trials (not 
necessarily in sequence) in which the submarine is detected, not sunk, does 
attack the FSE, and either sinks an FSE ship and is not killed or does not 
sink an FSE ship. The last trial is also preceeded by j - i - 1 trials in 
which the submarine is detected, not sunk, and does not attack the FSE. 
Denote the probability of each type trial, in the order listed, by Tl, T2, 


and T3. Then 


TG) G=— SRC) 
T2(1,K) = (PDFS(1,K) + PDFER) x (1 - PDFCE(I,K)) x (2 - FAIRSC) 
XJPSESE( I) xs (eee so IK): <P Sebo (ISK) 
and 
T301,K) = (PDE s(1,k) + PDRER) x (1 = PDRCE(1,K)) x (l= (2 =sEATRse) 


x PSFSE(I)) 


Therefore, 


-i-1 


Pe(C., D,) = [277] x (T2(1,k))* x (T3(1,K))9 2 x PTCT,K) 
bees at 


and equation (3.19) becomes 


: n+l 
a n RE Galicinoe tue 
E1(1,K) = ASALVO(I) yy ———————————————— 
=a i = (PSK 
1 = P9(I,K) 
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Penetration of the inner screen (ASW). Having penetrated the FSE, an 
attacking submarine must penetrate the inner screen prior to reaching the 
cargo ships. Figure 6 is the simplified tree of that part of Event 19 
(Submarine versus convoy) of SEALIFT from which the ASW associated probabil- 
ity of penetration is derived for this expected value analysis. This portion 
of the analysis is based upon an undated and informal modification of Event 
19 which, however, modifies only the ASW penetration portion. The assump- 
tions upon which the tree is constructed are the same as used for the FSE 
penetration analysis. That is, the ASW-detects-submarine probability is 
assumed not degraded by prior ASW ship losses and that if a submarine attacks 
the ASW, she is assumed to have sufficient torpedoes to do so. Because of 
the method used in deriving the ASW penetration probability for this analysis, 
parasitic torpedo usage could not be determined. It was observed that two 
mutually exclusive events prevailed, either the submarine penetrated the ASW 
or she was sunk. The method used for computing the probability of these 
events occurring was to first determine the probability the submarine was 
sunk. The desirtéd probability of penetration is found by complementation. 

The various paths to "submarine death" during the ASW penetration are 
stylized in Figure 7. The circled integers, (1), and the circled arabic A, 
(a), in Figure /7 represent branch points corresponding to those of the tree, 
Figure 6. The circled arabic X, (x), represents the event "submarine is 
killed." The symbol pe is the probability of transition from branch point 
r to branch point s, indexed by the integer i. From Figures 6 and 7 the 
probabilities for an (1,K) submarine are read off: 


PDETU(I,K) x VA(K) 


Pia 


PBB(I,K) x PAT(I,K) 


uU 
it 


A8 
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IS SUB DETECTED BY BOTTOM BOUNCE? 


YES NO 
PBB(I,k) | 
DOES SUB ATTACK UNALERTED INNER 
SCREEN? 
YES NO 
PSUBAT (I) 
IS A SCREEN DOES UNALERTED 
SHIP SUNK? INNER SCREEN 
wae ir DETECT SUB? 
PKOSC (1) YES NO 
PDETU(I,K) 
MR=) 


DOES HIGHLY ALERTED ® & 


INNER SCREEN DETECT 
SUB? 


YES NO 
PDETU(I,K) 


x VM(K) QTU(389) = 1? 


MR—=1 (3) ] ves 
S ( MR=-1 





DOES ASW SCREEN 


ATTACK SUB? 








YES NO 
PAT(1,K) Legend: 
BB=1 , 
DOES ALERTED SCREEN J: Submarine penetrates 
(8 ) DETECT SUBMARINE? ASW 
X: Submarine is killed 
= a INTEGER: Branch point - 
PDETU(I,K) 
x VA(K) 
MR 


FIGURE 6 


ARBORESCENCE ILLUSTRATING THE SUBMARINE-ASW SCREEN INTERACTION 
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IS SUB KILLED? 


YES NO 
PKBB(I,K) 


DOES SUB ATTACK THIS 
GROUP? 






YES 
PABB(I ,K) 






DOES SUB KILL A 
SCREEN SHIP? 


YES NO 
PKBBG (I) 








QTC(177) = 1? 


DOES THE SCREEN COUNTER ATTACK? 


YES NO 
PCC(1,BB) 


6 


Be = a7 


YES NO 


FIGURE 6 (continued) 
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DOES SCREEN KILL SUB? 


YES NO 
PKINS(I) 


DOES SUB ATTACK SCREEN? 










DOES SUB SINK A 
SCREEN SHIP? 


les NO 
PKSC (I) 


QTC(177) = 17 





DOES SCREEN COUNTER ATTACK? 


YES NO 


PSCA(I,MR) 


MR=MR + 1 


YES NO 


FIGURE 6 (continued) 
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DOES ALERTED SCREEN 
DETECT SUB? 






YES 
PDETU(I,K) 
x VA(K) 


FIGURE 6 (continued) 


4S 


Be 
- 








Pax (4 


Pax (4) 


FIGURE 7 


STYLIZED PATHS TO SUBMARINE DEATH DURING THE SUBMARINE~-ASW SCREEN INTERACTION 
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Poy = PKBB(I,K) 
Dae = (1 - PKBB(I,K)) x PCC(I,i) x (1 = PABB(I,K) + PABB(I,K) x 
(PKBBG(I) + (1 - PKBBG(I)) x (QTC(177)))), isl, ..., 3 
5 2 i i oe 
Pa) = l- Pop = Poy? i=l, e009 3 
4 
Pee PKBB(I,K) 
pa, 7 (1 = PKINS(I)) x PSCA(I,i) x (1 - PAA(I) x (PKSC(I) 
Pel, @ PKoeel jx ( PlesOreC L(y) Gas eed 
Pay = PKINS(I) 
rei SS PBB(I,K)) x PSUBAT(I1) x (1 - PKOSC(I)) x (1 - QTU(389)) 
and 
Pa, = PBB(I,K) x (1 = PAT(I,K)) x PDETU(I,K) x VA(K) 


+ (1 - PBB(I,K)) x (PSUBAT(I) x (PKOSC(I) x PDETU(I,K) x VM(K) 


+ (1 = PKOSC(I)) x QTU(389)) + (1 - PSUBAT(I)) x PDETU(I,K))) 


The probability that the submarine penetrates the inner screen, P3(I,K), can 
then be found as 
P3(1,K) = 1 = Pr (Submarine is sunk attempting ASW penetration) 


1 - Pr (a transit from the first (4) to (x) by all 


possible paths) x Pr (a transit from (A) to first 


(4) by all possible paths) + Pr (a transit from 


(A) to (x) by all possible paths not passing 


through any (4) ) 


47 





Each of these probabilities can be read from Figure 6 to be the following: 


Pr (a transit from the first (4) to (x) by all possible paths) 


Jee, eae oe 
Pax * Pag ‘Pax * Paq ‘Pax * Pag Pax 


)) 


Pr (a transit from (A) to. forse (4) by all possible paths) 


and 


it toe. Dens 
Par Pig t Pag!Pay Pig * Pag'Par Pia * Pag'Pa1 Pia * Pag Pai Pra 


> 


Pr (a transit from (A) to «) by all possible paths not passing through 


any (4) ) 


pi. (Pay + PanlPoy + ps (payne p> p 
ap "sx "ep ''ex ” Fes‘ ’ex ” “sas ex 


))) 


Convoy attack phase. It is particularly in the convoy attack that the 


flexibility of the SEALIFT game complicates the expected value analysis. 


Possible realizations of this phase depend heavily upon the chosen war 


option as listed in SEALIFT I 2 | and upon the various torpedo salvo sizes 


chosen for each submarine type. In view of these considerations and the 


problem which prompted this expected value analysis, a particular character- 


ization of the convoy attack phase has been chosen. It is assumed that 


1. 


one type of torpedo will be used against both cargo and 
screen ships, 
the salvo size against cargo ships is CSALVO(I) and that 


CSALVO(1) 2 


+ the initial torpedo supply, 
the salvo size against screen ships is ASALVO(I) and that 
ASALVO(I) 21/6 x CSALVO(I), 


in firing one CSALVO against cargo ships, the combination 


of CSALVO(I) and the probability of kill per torpedo 
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essentially eliminates the probability of sinking no cargo 
ships, and 

5. having fired two salvos at cargo ships, the submarine's 
torpedo supply is exhausted. 

Based on these assumptions, Figure 8 is the simplified tree of that 
part of Event 19 (Submarine versus convoy) of SEALIFT from which the convoy 
altack associated probabilities are derived from the first Markov chain in 
this analysis. For an (1,K) submarine sonar conditions the conditional 


probabilities of interest are 


P41(1,K) = Pr (a submarine shoots one salvo at cargo ships in the 
first Markov chain) 
= PDSAU(I,K) 
P42(1,K) = Pr (a submarine shoots two salvos at cargo ships in the 
first Markov ey 
= 1 ~ PDSAU(I,K) 
P5(I,K) = Pr (a submarine escapes/ she shoots one salvo at cargo 
ships in the first Markov chain) 
= 1 - PSSAU(I,K) 
and 
P6(I,K) = Pr (a submarine escapes/ she shoots two salvos at cargo 


ships in the first Markov chain) 


1 - PDSAU(I,K) x PSSAU(I,K) 

If the expected number of torpedoes fired by an (1,K) submarine at 
cargo ships in the first Markov chain given the submarine fires one and two 
Salvos are £4(1,K) and E5(1,K) respectively, then £4 and E5 are roughly 
approximated by: 


E4(I,K) 2 CSALVO(I) 
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SUB SHOOTS AT CARGO SHIPS 


’ 


IS SUB DETECTED BY SAU? 


NES NO 


PDSAU (I,K) G 





IS SUB SUNK? | SUB SHOOTS AT CARGO SHIPS 
YES No 
PSSAU (I,K) IS SUB DETECTED BY SAU? 
@ NO YES 
PDSAU(I,K) 
DOES SUBMARINE 
COUNTER ATTACK? SUB ESCAPES SUB SUNK? 
DECISION PSSAUS(I,K) 






ATTSAU(I) = 1 







SUB SINK A SAU 
SHIP? 


YES NO 
PSAUS(1,K) 


SUB ESCAPES 







SUB REATTACK CONVOY? 


NO YES 
NUMBER OF ATTEMPTS IS 
A PRE-SET DECISION 
GIVEN BY RASW(1I) 


SUB ESCAPES NOT CONSIDERED 


FIGURE 8 


ARBORESCENCE ILLUSTRATING THE CARGO SHIP ATTACK PHASE 
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and 
E5(1,K) = TORPS(I,ASW) - CSALVO(I) - E1(I,K) 
The tree aSsociated with the cargo attack phase for those submarines ~ 
which successfully enter the second Markov chain is obtained from Figure 9 


by changing the entire right branch at the branch point K into the event 


"Submarine Escapes." The conditional probabilities of interest are then 


P4(1,K) = Pr (a submarine shoots her remaining torpedoes at cargo 
ships) 
= 1 
Pll(I,K) = Pr (a submarine escapes given she shoots) 


Poth.) 
The expected number of torpedoes a submarine fires at cargo ships in 
the second Markov chain given she fires, E6(1,K), is roughly approximated by 


E6(I,K) 3 TORPS(I,ASW) = CSALVO(I) - 2 x El1(I,K) 





The transition probabilities. The calculation of the Markov chain 
transition probabilities needed for this expected value analysis can now be 
obtained by combining the conditional probabilities derived in the previous 


sections. fhese transition probabilities are: 


Paae > 6 lee, P1(I,K) 

Po = Gy = (1 - P1(I,K)) x P9(I,K,n) 

P3 = (1 - Pl1(I,K)) x P7(I,K,n) x P3(I,K) x P41(I,K) x P5(I,K) 

Py, = (1 = P1(1,K)) x PY(I,K,n) x PA(1,K) x (1 = P5(1,K)) 

Po = (T= P11 jK))-« P71 Kf) x P30, K) x P4201 Ke x Peek) 

Pe = (1 — PL(I,K)) x PT(1,K,n) x P3(1,K) x P42(1,K) x (1 = P6(1,k)) 
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Pp = G6 = (l= Pik) x PP tl,K,n) x (1 = P3C1,K)) 

Pg = Ag = (1 - PL(I,K)) x PB(I,K,n) 

q, = (1 - P1(1,K)) x P7(1,K,n) x P3(1,K) x P4(1,K) x P5(1,K) 

1 = (1 - Pl(I,K)) x P7(1,K,n) x P3(1,K) x P4(1,K) x (1 = P5(I,K)) 


3.2.3 Expected Number of Torpedoes Fired at Cargo Ships 
In order to determine the expected number of torpedoes fired at cargo 


ships for any given cycle, define the following to be 


Pine) 2 event an (1,K) submarine fires u salvos at cargo ships 
in the v th Markov chain, u, v =1, 2 

NS(w,1,K) : the number of (I,K) submarines which arrive on station 
for the w th cycle, wi= 1, 2, 3,200 | 

NT(w,1,K) the total number of torpedoes fired at cargo ships during 


the w th cycle by (1,K) submarines, w = 2, 3, eee 
The following equation is then used 
(3.20) | 
E(NT(w,1,K)) = E(NS(w,1I,K)) x [Pr(F(1,1,1,K)) x E4(1,K) 
+ Pr(F(2,1,1,K)) x E5(1,K) 
+ Pr(F(1,2,1,K)) x E6(1,K) 
For the first cycle, distinction is made for the initially on-station, 
initially not-on-station contingency. For each submarine type and for each 


sonar condition the expected number of torpedoes fired in each contingency 


is computed as 


E(NT, (1,1,K)) E (number of torpedoes fired by (I,K) submarine initially 


on station during the first cycle) 


aye 





E2(1,K) x | PA2(1,K) x E4(1,K) 
+ PA3(I,K) x E5(1I,K) 
+ PAS(I,K) x E6(1,K) 


and 


E(NT.(1,1,K)) E (number of torpedoes fired by (1,K) submarines 


initially not-on-station during the first cycle) 


E3(1,K) x POC(I,K) x9 PAT(I,K) x E4(1,K) 
+ PA8(I,K) x E5(1,K) 


+ PALO(I,K) x E6(I,K) 


The expected number of survivors from the first cycle, i.e€., those who safely 


reach home port, EC(1,I1,K), is 
EC(1,1,K) = (E2(1,K) x PA4(I,K) + E3(1,K) x PO(I,K) x PA9(I,K)) x PO(I,K) 


According to the rules of SEALIFT, the survivors are re-apportioned among 
the sonar boxes to maintain balance. The expected number of re~apportioned 


survivors of the (I,K) classification, EC(2,1,K), is 
EC(2,1,K) = EC(1,1,K) x > PIT BOX( J) 


where the sum is over those sonar boxes of like sonar conditions. To deter- 


mine E(NT(2,1,K)), equation (3.20) is applied again, and becomes 


EXNTC2515K)), = ae ne E(NT.(1,1,K)) 


The expected number of survivors from the second cycle, EC(3,1,K), is 
EC(3,1,K) = EC(2,1,K) x PA9(I,K) x PO(I,K) x PO(I,K) 


For the w th cycle, w =3, the expected number of re-apportioned survivors 
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from the previous cycle, EC(2w-2,1,K), serve as the submarine input, and 


the following iteration formulae are found, 


EC(2w-3,1,K) x EC(2w-4,1,K) 


EC (2w-2,1,K) EC (2w-5,1,K) 


and 


E(NT(w-1,1,K)) 
EC(2w-4,1,K) 


E(NT(w,1I,K)) EC(2w-2,1,K) x 


The expected number of survivors from the w th cycle (not yet re-apportioned 


among the sonar boxes), EC(2w-1,1,K) is 


; EC(2w-3,1,K) 
EC(2w-1,1,k) = EC(2w-2,1,K) x EC 4 TK) 


An estimate of the number of cycles for the war is given by the minimum 
number of cycles a live submarine, not initially on station, might expect to 
achieve. This number is a function of submarine type and is given by the 


integer L(I) 


ean TWAR 
. DSUBSTA 
ALTIM(I) + DTSPRT(I) + Sash 


An estimate for total expected number of torpedoes fired at cargo 
ships by (I,K) submarines for the war, E(NT(I,K)), is then 


L(1) 


E(NT(I,K)) = E(NT, (1,1,K)) + E(NT.(1,1,K)) + E(NT(w,1,K) 
w=2 . 


Multiplying this number by the single shot kill probability, PSCS(I,K), pro- 
vides an estimate for the number of cargo ships sunk by (I,K) submarines, 


E(CS(I,K)), 


E(ES(1{K)) => ECNE(1 5K) ) xsPSES(1,K} 
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Finally, summing E(CS(1,K)) over the indicies I and K provides an estimate 


for the total expected number of cargo ships sunk during the entire war, E(CS) 


3 
E( ES) = > \ E(ESGIGK)). 


[zl Ke=l 


ala 





4. Conclusions and Recommendations 

As a model purporting to represent a probabilistic real world situation, 
the SEALIFT game, the analysis presented in this thesis remains largely un- 
tested. The opportunity and data base for testing the model's reliability, 
accuracy and sensitivity to the proposed approximations are, however, readily 
available. If the SEALIFT game is contemplated as an analytical tool in the 
CYCLOPS III study, the accuracy of this expected value model should be ex- 
plored. An acceptably reliable programmed version of this model should 
effect appreciable savings in computer time. 

For her continued patience, encouragement and good humor, I acknowledge 


the support of my wife in the preparation of this thesis. 
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APPENDIX I 
DETAILED DESCRIPTION OF THE SEALIFT GAME 
The following description is quoted from SEALIFT I {| 


II. GENERAL DESCRIPTION 


Convoys 


In this game, friendly forces have two ports to transit between home 
port and delivery port. During this transit the convoy may cross up to two 
"constant" barriers (oneway transit). The constant barrier in this simulation 
can sink ships that cross it but cannot itself be weakened by these sinkings 
or by the actions of the ships crossing it. Mine fields, air barriers, etc., 
can be described as "constant" barriers. The "variable" barrier, used against 
submarines, can sink ships that traverse it but can be weakened by having the 
ships that comprise it sunk. Barriers of this type might be submarine bar- 
riers, surface ship barriers, etc. The "variable" type of barrier will be 
described in greater detail later in this section. It should be pointed out 
that while each barrier that could be encountered by a convoy or submarine 
will be described herein, any one or all of these barriers may be removed for 
a desired play simply by setting to zero its probability of detecting a pro- 
spective "transitor"... 


At certain times, which are inputs, a convoy is formed in the home port, 
The convoy is made up of three types of ships: cargo ships, ASW escort ships, 
and forward screen ships*. Each convoy when it forms draws its required 
complement of ships from a "port pool." This pool is stocked by ships from 
returning convoys and by the ship building rate. If the required ships are 
available, the convoy sails toward its delivery port: if not, the game is 
stopped, and a "printout" is made of the game results to date including the 
event that caused the stoppage. 


At one hundred miles from its home port, all ships of the convoy except 
the forward screen escorts pass through a constant barrier (barrier 1C€). The 
losses, if any, are recorded and subtracted from the convoy. 


After having transited the first constant barrier, the convoy becomes 
liable to an air attack. This attack is played as the probability of a con- 
voy ship being sunk per day at sea by air attack. The convoy then enters 
the submarine patrol area and may encounter a submarine while in this area. 
The submarine/convoy encounter will be described in greater detail later. 


Ships lost because of any one of the aforementioned reasons are noted 


by type of ship sunk, method of sinking, and time of loss, and are subtracted 
from the convoy to be printed out at the close of the simulation. If these 


* Variations in inputs permit the forward screen to be composed alternately 
of aircraft, sonobuoys, submarines, etc. 
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losses include ASW escorts or forward screen escorts, the convoy's ability 
to detect** enemy submarines will be diminished. Each time ASW escorts 
and forward screen escorts are sunk, new lower detection probabilities for 
the various parts of the convoy are calculated. 


After passing through the submarine patrol area, the convoy becomes 
liable to another air attack, then passes through a second constant barrier 
one hundred miles from its delivery port. After reaching the delivery port, 
the convoy unloads and stays in port a prescribed number of days. During 
this period the cargo ships are again liable to air attack. Of those ships 
sunk in an air attack, one-half are assumed to be sunk before unloading and 
one-half after unloading. [The number of cargo ships delivered are now 
tabulated and the convoy is ready to proceed home. 


Proceeding home, the convoy follows the same route in reverse order, 
passing through the same barriers, the same submarine patrol area, and the 
Same series of air attacks. When it reaches its home port, the number of 
ships remaining in the convoy are placed in a "ship pool" by ship type. 
These ship pools will then be used to make up future convoys. 


Submarines 


This model is designed to handle three types of submarines. All of the 
parameters which describe the submarine types are inputs to the game. At 
the start of the war, the opposition submarines, chosen by type of submarine 
as well as by percentage of each type, are placed in two predetermined "game 
positions," selected according to the option of the user. The two game po- 
Sitions are "on station" and "in port." For those submarines that start the 
game in port, a given number (input) will leave at a given time interval 
(input) until all have left for their stations. 


While a submarine is in port it may be susceptible to strikes at source. 
In addition, the submarine building rate creates new submarines as the war 
progresses. As each new submarine or "turnaround" submarine leaves port, 
it is given a full load of anti-ASW and anti-shipping torpedoes, and a time 
at which it must depart from station and return home. Both of these inputs 
are a function of submarine type. Submarines that are on station at the 
Start of the war are given a full load of both types of torpedoes, but their 
time remaining on station is determined by a uniform random distribution, 
The enemy submarine patrol areas are considered to be a constant distance 
from the submarine port. 


While all submarines are considered to transit the same average distance 
to their patrol area, they can appear in one of four areas along the convoy 
route, each area involving either good or bad sonar conditions. Each area 
is looked upon as a "sonar box” and is given a certain percentage of each 
type of submarine (input). A convoy must pass through the sonar boxes in a 


** The terms "detect" or "detection" as used throughout this research contri- 
bution refer to the process of sensing a signal and classifying it as due 
to a target. 
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fixed sequence, viz.e, Box 1, first, Box 2, second, etc. Either good or 

bad sonar conditions may be assigned to each box. If a given submarine is 
put into a certain box, the submarine derives both its detection and sinking 
probabilities with respect to the convoy, not only from its type, but also 
from the sonar conditions assigned to the box. Thus for each decision point 
in the submarine/convoy engagement, there is in effect a 2 x 3 probability 
matrix, the rows of which correspond to the two possible sonar conditions, 
and the columns to the three possible types of submarines. 


When the submarine leaves port, it may be required to pass through as 
many as eight barriers, four constant and four variable, before arriving in 
the submarine patrol area. Each barrier can be used or not used as desired 
(input). The constant barriers are of the same type as those described for 
the convoy, except that in addition to having a probability of being sunk 
as it transits the barrier, the submarine can also incur a time delay as it 
goes through or around the barrier. Any time delay so incurred will of 
course cut down the remaining time the submarine can spend in its patrol 
area. 


When the submarine crosses a variable barrier, it not only has the chance 
of being sunk itself, but also has the chance of sinking a ship in the barrier. 
If the submarine does sink a ship in the barrier, the probability of detection 
for this barrier is assigned a new lower value. The barrier ship sunk is 
replaced at some later time in the game (an input), and the detection prob- 
ability of the barrier increased. In this way the variable barrier not only 
destroys some of the submarines that must transit it, but is itself sub-= 
jected to possible reductions in effectiveness. Of course, submarines which 
are detected by the barrier but not sunk suffer time delays which reduce 
the time that can be spent in the patrol area. 


As long as the submarine is at sea, it becomes susceptible to being found 
and engaged by a HUK group. There are two different probabilities per day at 
sea that a submarine will encounter a HUK group; one applies when the sub- 
marine is in transit to or from its patrol area, and the other when the sub- 
marine is in its patrol area waiting for a convoy. Should the submarine be 
found by a HUK group, the HUK group will attempt to sink the submarine. In 
turn, the submarine has a chance of sinking one of the HUK ships. 


The time of arrival "on station" depends upon the speed of the submarine 
(which in turn depends upon submarine type) and upon the time delays the 
Submarine undergoes while crossing barriers and fighting HUK groups. fhe 
submarine remains on station or in its patrol area (sonar box) until a con- 
voy crosses the "subline." The length of the "subline" or what could better 
be called the submarine patrol area, is an input. 


It is assumed that a convoy crosses the center of the submarine patrol 
line. When this happens, the first submarine is drawn from Box 1 and 
placed, using a uniform random distribution, on the subline. Its position 
on this line, relative to the center of the line, is compared to an athwart- 
ship detection range for this type of submarine and for its sonar box. If 
it 1s within the detection range, the submarine detects the convoy, other- 
wise it does not. If it detects the convoy, it will attempt to "close" on 
this convoy; if it does not detect the convoy, it will be returned to its 
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sonar box (patrol area). 


The first three boxes are played in this manner: the fourth box is 
played slightly differently. As the convoy approaches the fourth box 
(fourth patrol area), there is given a probability per month of the war 
that submarines in this patrol area will contact the convoy. If contact 
is made, this box will be played as the other three. If, however, this 
contact is not made, this box will not be played for this convoy. This 
fourth box might then be used for a submarine wolf pack; for, although a 
wolf pack has a low probability of intercepting a convoy, given this in- 
terception it has a high probability of sinking ships in the convoy. 


In general, if a submarine has detected a convoy, it will next try 
to penetrate the convoy defenses and sink cargo ships. (However, an option 
is permitted, if desired, by which submarines will try to sink ASW escorts 
in an anti-ASW escort campaign.) The event that plays the submarine en- 
gaging the convoy is Event Subroutine No. 19, "Submarine vs. Convoy." In 
this event, the submarine must first penetrate the forward screen (or out~ 
ward screen of the convoy). Given a successful penetration of this screen, 
it must then penetrate the ASW screen (or innermost convoy screen). When 
it has successfully penetrated both of these screens, it may then attack 
cargo ships. When the submarine has sunk one or more ships, either cargo 
or screen ships, the surviving screen ships form surface attack units 
(SAU's) which attempt to locate the submarine and attack it. Anti-shipping 
and anti-ASW torpedoes are used during this attack, and whenever the allow- 
able minimum number of torpedoes of the type needed to continue the attack 
has been reached, the submarine is required to break off the attack and 
start its transit for home. ... 


After a submarine has either reached its minimum number of torpedoes 
or has used its allowable time, whichever occurs first, it must return to 
its home port. Qn the way home, the submarine must traverse, in reverse 
Order, the same barriers that it crossed on its way out to station. When 
the submarine arrives at its home port, it will stay there for a certain 
number of days (input) and then be resupplied and again sail for its patrol 
areas 
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